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(Abstract) This paper proposesan improved load flow method based on perturbation theory which is simple, efficient and reliable. 
The performance of the proposed method has been tested on well-conditioned IEEE 5-bus, 14-bus, 30 bus and also tested on 
11 -bus ill conditioned test systems. The proposed method shows its efficacy and robustness with minimum number of iterations 
along with high accuracy on conventional load flow methods like G-S (Gauss-Seidel) method. 
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1. INTRODUCTION 

Load flow calculations are imperative for operational planning, 
control, system stability and reliability. The load flow study 
[l]-[4] is an important tool involving numerical analysis 
applied to a power system. Study of load-flow studies is very 
important as present electrical power system which is 
continuously expanding in size and growing in complexity. 
This study ultimately gives important conclusions while to 
determine the steady state operation of an Electric power 
system. 

The principal information obtained from the power flow 
study [5]-[6] is real & reactive power flowing in each line, 
voltage magnitude, phase angle at each bus. The development 
of digital solution can be found in [7]-[9]. There are various 
methods are available and some modifications on the existing 
methods are shown in [10]-[1 1]. In this paper a simple, efficient 
and reliable method based on Power Perturbation has been 
described which involve lower storage requirement along with 
advantages like convergence and accuracy of conventional 
techniques. As a result it can also be used in real time 
applications. 

2. LOAD FLOW PROBLEM 

It is convenient to represent power system networks using the 
single line diagram, which can be thought of as the circuit 
diagram of the per -phase equivalent. In general, representation 
of a single line diagram is used for load flow studies. 

Figure 1 shows the schematic diagram of the basic structure 
of a typical bus system. It consists of n number of lines with 
lumped parameters like line voltages. Vi, V 2 , ....V„ and line 
admittances like Y n , Y i2 , Y in respectively and Power node 



Equations and subsequent Mathematical model for the said 
problem can be enumerated as follows: 



Figure 1 : Basic schematic diagram of a typical bus system, 
i) Current injections at any node n can be written as, 

h = y m v t + y n (y, -v t ) + y a <y, -v 2 ) + ... + yjv, - v„) 
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ii) The real and reactive power at bus i is 

p i +jQ i =y i i- o) 
it=m-jQiVv;] (4) 

Power flow equations formulated in polar form. For the 
system in Fig.l, Eqn.2 can be written in terms of bus 
admittance matrix as: 
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Expressing in polar form: 
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Substituting for //from Eqn.6 in Eqn. 4 
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Separating the real and imaginary parts, 

P i=f.\Vi\VM c ™%- S i+ S ^ (8) 
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Expanding Eqns. 8 & 9 by Taylor's series results to 
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3. GAUSS SEIDEL METHOD 



The Gauss-Seidel method [12]-[14] is an iterative process 
which starts by assigning estimated values to the unknown bus 
voltages. Using the estimated bus voltages and the specified 
real and imaginary power values, a new value for each bus 
voltage is calculated at the end of each iteration. The process is 
repeated until the difference between each bus voltage and its 
corresponding value in two successive iterations is less than a 
predefined tolerance value. 

For a total of N buses, the calculated voltage at any bus K is 
expressed as a function of the real and reactive power delivered 
to a bus from generators or supplied to the load connected to the 
bus, the estimated or previously calculated voltages at the other 
buses, and the self- and mutual admittances of the nodes. 
Equation (5) is an algebraic non linear equation which must be 
solved by iterative techniques and V ; can be obtained by solving 
iteratively. 
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Where yy is the actual admittance in p.u. the specified real & 
reactive power (P isp &Q iS p) are in p.u. by assuming current 
entering bus I was positive. Thus for Generator buses these 
injected real and reactive powers have positive values. For 
Load buses the real and reactive powers flow away from the bus, 



-2- 



P isp and Q isp have negative values. By solving Equation (5) 
the P; and Q ; are 
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The power flow equation is usually expressed in terms of the 
elements of the bus admittance matrix, Y bus _ shown by upper 
case letters, are Yy = -yy, and the diagonal elements are Y ;i = X 
yy. Therefore Equation (6) can be written as 

pr - mr 

y*(k) 



+ 



Yy..v { 

L-t 'J 7 



V k+1 =- 



,J*l 



(14) 



pk+l = R (y* W { y* W ^ yW j }j ^ . (15) 



7=0 



7=1 



or = -l (y? k) w; ik) i y v - x y, - v™ (i6) 

7=0 7=1 

Iterative Stages for the above method is as follows: 

Slack bus: both components of the voltage are 
specified. 2(n-l) equations to be solved iteratively. 

• Flat voltage start: initial voltage of 1.0+j0 for 
unknown voltages. 

• PQ buses: P isp and Q isp are known, with flat voltage 
start, Equation (9) is solved for real and imaginary 
components of Voltage. 

PV buses: P isp and [Vi] are known. Equation (11) is 
solved for Q k+1 which is then substituted in Equation 
(9) to solve for V k+1 

However, since [VJ is specified, only the imaginary part of 
V; k+1 is retained, and its real part is selected in order to satisfy 

(e 1 (k+1) ) 2 + (f 1 (k+1) ) 2 =IV,l 2 
(17) 
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Acceleration factor: It is the rate of convergence is increased by 
applying an acceleration factor to the approximate solution 
obtained from each iteration. 

V* +1) = V, (k) + a (V Ical (k) _ V/ k) ) (19) 
Iteration is continued until 

| ei (k+1) . ei (k) l <8, |fi (k+1) . f, (k) l <8 (20) 
Once a solution is converged, the net real and reactive powers 
at the slack bus are computed from Equations (10) & (1 l).In the 
original algorithm of Gauss-Seidel method the voltage at each 
bus at any iteration is calculated by using the voltages of 
previous buses calculated at the same iteration and the voltages 
of the next buses calculated at the previous iteration i.e. where 
the difference of corrected voltage precision index and previous 
voltage is within the acceptable limit. The entire process is 
carried out again and again until the amount of correction in 



voltage at every bus is less than some predetermined precision 
index. 

It has been professed that the G-S method for solving large 
power flow problems takes excessive number of iterations and 
more computational time. However the proposed Power 
Perturbation method shows that the total execution time is 
considerably reduced. 

4. POWER PERTURBATION METHOD 

Very often, a mathematical problem cannot be solved exactly 
or, if the exact solution is available, it exhibits such an intricate 
dependency in the parameters that it is hard to use as such. 
Perturbation theory gives a systematic answer to this question. 
The proposed methodology can be outlined by the flowchart as 
shown in Figure 2 
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Perturbation theory [15]-[18] comprises mathematical 
methods that are used to find an approximate solution to a 
problem which cannot be solved exactly, by starting from the 
exact solution of a related problem. Also this theory leads to an 
expression for the desired solution in terms of a power series in 
some "small" parameter that quantifies the deviation from the 
exactly solvable problem. The leading term in this power series 
is the solution of the exactly solvable problem, while further 
terms describe the deviation in the solution, due to the deviation 
from the initial problem. Formally, we have for the 
approximation to the full solution A, a series in the small 
parameter (here called e), like the following: 
A=£°A +e 1 A 1 +e 2 A 2 + 



A would be the known solution to the exactly solvable 
initial problem and A ly 4 2 ,... represent the higher order terms 
which may be found iteratively by some systematic procedure. 
For small e these higher order terms in the series become 
successively smaller. An approximate "perturbation solution" 
is obtained by truncating the series, usually by keeping only the 
first two terms, the initial solution and the "first order" 
perturbation correction: A= A +eAi 

5. PERFORMANCE ANALYSIS 

The performance of the P-P Method was tested on IEEE 5-bus, 
14-busand 30-bus as well as on ill conditioned bus system [19] 
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with a convergence accuracy of 10 on a MVA base of 100 or 
-i 

equivalently 10 MVA for both power residuals AP and AQ. It 
requires the same storage size as the Gauss-Seidel method, but 
less computation time than the conventional Gauss-Seidel 
program. It is attractive for accurate or approximate results and 
contingency calculations for IEEE test systems. The proposed 
solution algorithm can be implemented directly and efficiently 
into microcomputers with limited storage capacities. To 
examine the effectiveness, test systems including IEEE 5-bus, 
14-bus, 30-bus and 11-bus ill conditioned bus systems were 
studied by the proposed method, and the convergence 
characteristics of these bus systems are compared with the 
Gauss-Seidel method. 

5.1 Number of Iteration 

The Gauss-Seidel method of solving large power flow 
problems it has been seen that an excessive number of iterations 
are required for different IEEE test systems compared to the 
proposed Power-Perturbation method before the voltage 
corrections are within an acceptable precision index. The test 
results are given in Table 1 and the corresponding graph is 
plotted in Figure 3. The power flow solution based on power 
perturbation method shows its robustness and assured 
convergence for all types of test systems with minimum 
number of iterations and also it is very important that P-P 
Method converged for an ill conditioned system (1 1-bus) where 
conventional G-S Method diverged. 



Figure 2. Procedure of Power perturbation method 
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5.3 Accuracy 



Figure 3.Comparison of G-S & P-P Method (Iterations for different 
IEEE bus systems) 

5.2 Computation Time 

The algorithm has been developed and run on Intel core 2 duo 
processor, 1.86 GHz, 2 GB RAM, 250 GB of Hard-disc laptop 
using MATLAB 7. The results of time-comparison are given in 
Table 2. In conventional Gauss-Seidel based power flow 
algorithm, Jacobian matrix is required. 

As memory requirement of matrixes is more, CPU time per 
generation for conventional method is also more. The power 
perturbation technique shows better 'time per iteration' 
compared to the conventional Gauss-Seidel method. 

The convergence time for G-S and proposed P-P method are 
shown in Table 1, where it is clearly depicted that G-S Method 
fails to converge for ill conditioned system where as P-P 
method gives the solution. 

The same is shown in the Figure 4 in a graphical form which 
reveals that time required for Power Perturbation method is 
smaller than the conventional G-S method. 
Table 1. Comparison of Computation Time in Second 
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Both G-S & P-P methods were tested by varying accuracy from 
10" 1 to 10" 3 the and was found that the Gauss- Seidel method 
will have maximum number of iterations compared to the 
proposed Power-Perturbation method as shown in Table 2. 
Table 2. Comparison of Accuracy and Iterations 
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Figure 4 Comparison of G-S & P-P Method (Convergence time in 
seconds for different IEEE bus systems) 



6. CONCLUSION 

Power Perturbation method gives faster, smarter and more 
accurate results pertaining to the conventional Gauss-Seidel 
technique of load flow calculation.Test results of the proposed 
method is performed on the well-conditioned IEEE 5-bus, 
14-bus and 30-bus system as well as on ill conditioned 11-bus 
system. 

The numerical results show that the P-P method is indeed 
feasible, effective and is superior to the G-S method. Further 
research work can be done on Power-Perturbation method for 
more complex bus system as well as ill-conditioned systems to 
verify its ability of convergence, robustness and the pace at 
which it is providing the solution. 
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